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Abstract. Observational and theoretical investigations provide evidence for non-uniform spot and magnetic flux distributions 
on rapidly rotating stars, which have a sign ificant imp act on their angular momentum loss rate through magnetised winds. 
Supplementing the formalism of MacGregor & Brenner 1 1991) with a latitude-dependent magnetised wind model, we analyse 
the effect of analytically prescribed surface distributions of open magnetic flux with different shapes and degrees of non- 
uniformity on the rotational evolution of a solar-like star. The angular momentum redistribution inside the star is treated in 
a qualitative way, assuming an angular momentum transfer between the rigidly-rotating radiative and convective zones on a 
constant coupling timescale of 15Myr; for the sake of simplicity we disregard interactions with circumstellar disks. We find 
that non-uniform flux distributions entail rotational histories which differ significantly from those of classical approaches, 
with differences cumulating up to 200% during the main sequence phase. Their impact is able to mimic deviations of the 
dynamo efficiency from linearity of up to 40% and nominal dynamo saturation limits at about 35 times the solar rotation rate. 
Concentrations of open magnetic flux at high latitudes thus assist in the formation of very rapidly rotating stars in young open 
clusters, and ease the necessity for a dynamo saturation at small rotation rates. However, since our results show that even minor 
amounts of open flux at intermediate latitudes, as observed with Zeeman-Doppler imaging techniques, are sufficient to moderate 
this reduction of the AM loss rate, we suggest that non-uniform flux distributions are a complementary rather than an alternative 
explanation for very rapid stellar rotation. 
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1. Introduction 

In the presence of open magnetic fields the angular momen 
turn (AM) loss of a star through winds and outflows is signifi 
cantly enhanced since the tension of bent field line s effectivel 
prolo ngs the lever arm of the associated torque dSchatzma 
1962). Whereas during the pre-main sequence (PMS) phase 
the rotational evolution of a star is dominated by its chang- 
ing stellar structure and magnetic interaction with a circum- 
stellar accretion disk, its rotation during the main sequence 
(MS) phase is mainly determined through braking by mag- 
netised winds jBelcher & MacGregoJl97a) . Theoretical stud- 
ies of magnetised win ds go back to class ical app roaches of 
IWeber & Davisl i 19671 hereafter WD) and lMestell ( 1 19681) . At 
high rotation rates the magnetic field adds considerably to 
the acceleration of the outflow through magneto-centrifugal 
driving (Michel 1969), which makes magnetised winds in- 
trinsically latitude-dependent. The total AM loss rate is con- 
seque ntly susceptible to variations of the surface magnetic 
field JSolanki et alJll997t lHolzwarthl|2005i) and atmo spheric 
field topology jMestel & Spruitlll987t lKawaleJll988l) . More 
recently, multi-dimensional MHD-simulations have been ac- 
complished to study structural and temporal wind properties 
like coronal mass ejections in more detail (e.g.. lSakurailll985t 



Kep pens & Goedbloedll999l l2000l) : the extensive computa- 
tional requirements render them however less attractive for 
studies concerning the rotational evolution of stars. 



1.1. Saturation limits 
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ISkumanicbl dl972l) analysed the rotation and chromospheric 
emission of cool stars in different evolutionary stages, and 
found that their rotation rate (as well as Call luminosity) is 
about proportional to the inverse square-root of their age, that 
is Q a r 1 ' 2 . Presu ming a const ant moment of inertia and a 
WD braking law, the Skumanich relation implies a linear rela- 
tionship between the rotation rate of a star and its characteris- 
tic magnetic field s trength, B o c Q" de with dynamo efficiency 
«de — 1- However, Saar (1991) found that the magnetic flux 
rather than the field strength is increasing linearly with the ro- 
tation rate, which is now commonly adopted (at slow rotation 
rates). The consequence of a continuously increasing stellar 
field strength is a very efficient magnetic braking of rapidly 
rotating stars. In fact, too efficient, since observations of young 
open clusters reveal significant nu mbers of stars with rotational 
velocities up to v sin i ~ 200 km/s (Stauff er et alll997l) . whose 
existence is, with initial rotation rates o f young T Tauri star s be- 
ing observationally well constrained ( B ouvier et alJ 1993), dif- 
ficult to explain in the framework of magnetic braking without 
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any mechanism which moderates the AM loss rate at higher ro- 
tation rates. Respective studies therefore frequently presume a 
saturation of the AM loss beyond a limiting rotation rate, whose 
value is very susceptible to model assumptions like the internal 
AM redistribution; where as some investigations require low (< 
20 Qp) saturation limits (jk eppens et aljjl995|: barnes & Sofia 
1 199ft Boil vier et alJll997l iKrishnam urfhi et allll997l), others 
argue for higher (> 40 Q ) va lues ( Soderbl om et alJ 1 19931 
ICollier Camero n & Jian kell994l) . 

The saturation of the total AM loss rate is often ascribed 
to the underlying dynamo mechanisms in the convective en- 
velope, because current dynamo theories anticipate a moder- 
ation and eventually saturation of the amplification process 
(rid e — > 0) when the back-reaction of strong magnetic fields 
supp resses plasma motions and inhibits its further increase 
(e.g. jRiidiger & Kichatino\ll993l) . However, the various theo- 
retical models are as yet unable to provide consistent or explicit 
values for the critical rotation rate or field strength at which this 
occurs. 

The concept of a saturation of magnetic activity is strength- 
ened by empirical activity-rotation relationships, which reveal 
(for rotation periods longer than a few days) close correla- 
tions between the rotation rate and the strength of activity 
proxies like magneticall y induced chromospheric and coronal 
emission JNoves et all 1 1984 IVilhulll984[ Mathiou dakis et alJ 



fl995tlStauffer et alJl997tipizzolato et alfcOOfl. In rapidly ro- 
tating stars several activity signatures are found to saturate: the 
chromospher ic UV emission below rotation periods P ~ 3 d 
JVilhd Il984h. and the EUV and (soft) X-ray emission for 
P < 2 d JStauffer et alJll997h . However, the variation of pho- 
tometric light curves, associated with the presence of dark 
spots in the stellar photosphere, is found to increase for even 
shorter rotation periods, down to P ~ 0.35 - 0.5 d, for which 
the c hromosp heric and coronal emission are already saturated 
dO'Dell et all Il995t iMessina et alJ l200lh: f or an antithetical 
point of view see Krishnamurthi et al J Jl998l) . 

The saturation of activity signatures is not unambigu- 
ously indicative of a saturation of the dynamo mechanisms 
in the convective envelope, since emission processes are li- 
able to further rotation-dependent effects like a reduction of 
the X-ray emitting volume through the centrifugal stripping 
of hot coronal loops, or the shift of coronal loop tempera- 
tures into different emission regimes as the effective gravi- 
tation and pressure scale hight change with the rotation rate 
JUnruh & J ardine 119971: Jardine & Unruhl ll999l) . The conjec- 
ture that changes of the atmospheric emission are not necessar- 
ily correlated with the (sub-)photospheric magnetic activity is 
supported, for example, by observations of the ultra-fast rota- 
tor VXR45a (V370 Vel, P = 0.223 d), whose X -ray emission 
is bel ow the typical level of X-ray saturated stars ( M arino et alJ 
120031) . whereas its brightness surface maps are stil l very simi- 
lar to those of more slowly r otating stars (Mars den et all2 004: 
IVrielmann & Hussainl20ol . 

The large range of rotation rates (~ 10-70 Q ) in which ob- 
served activity signatures are found to saturate raises the ques- 
tion whether these phenomena reflect the actual behaviour of 
the underlying dynamo processes, in particular beyond which 
critical rotation rate their efficiency breaks down. Since the 



range of observed saturation limits practically covers the one 
of suggested AM loss limits, a definite justification of the latter 
in terms of a dynamo saturation is rather precarious. 

1.2. Surface magnetic flux distributions 

Doppler imaging (DI) observations of rapidly rotating stars 
yield non-uniform surface brightness distributions, where, in 
contrast to the case of the Sun, dark spots are not only located in 
equatorial regions, but also at intermediate and polar latitudes 
JStrassmeieill2002( and references therein). Theoretical mod- 
els considering the formation of magnetic features at higher 
latitudes involve the pre-emptive poleward deflection of mag- 
netic flux inside the conve ction zone by the Coriolis force 
JSchiissler & Solanki Il992h . and/or its post-eru ptive pole- 
ward transport through meridional motions JSchriiver & Title! 
200 ll). Backed by these observational and theoretical results, 



Sola nki et all l l 1997b investigated the influence of a bi-modal 



magnetic field distribution on the rotational evolution of cool 
stars. They found that a concentration of magnetic flux at very 
high latitudes reduces the total AM loss rate as efficiently as a 
dynamo saturation limit at ~ 20 Q . Based on their findings, 
they question the concept of a dynamo saturation at low rota- 
tion rates and argue instead for a saturation above 50 Q ; a sim- 
ila r though more qualitative argument has also been discussed 
bv lBuzasillll997l) . 

The work of Solan ki et all was focused on flux concentra- 
tions around the pole, assuming that the observed brightness 
distributions indicate likely locations of open field lines along 
which a stellar wind can escape. But the bare existence of 
starspots does not a priori imply information about the associ- 
ated magnetic field topology, that is neither about the existence 
nor about the amount of open magnetic flux. Zeeman-Doppler 
imaging (ZDI) observations in contrast directly confirm the 
magnetic origin of the dark features and enable a determination 
of the magnetic field at the stellar surface dDonati et al.lll997T) . 
Such field distributions s erve as boundary c onditions for field 
extrapolation techniques jAltschuler & Newkkklll969l which 
reveal large-scale magnetic field topologies and consequently 
the d istributio n of clo sed and open magnetic field lines (e.g., 
Llardine et alll2002albh . In the case of the rapidly rotating star 
LQ Hya (P = 1.6 d) the latitudinal distributions of open flux 
show tha t large amounts are located at intermediate and high 
latitudes (Mclvoret al. | |2004 . A similar result is found in the 
case of A B Dor (P _ 0.51 d) for observations between 1995 
and2003 dDonati & Collier Cameronll 997llDonati et al Jl 99^ 
2003): Figure[0shows the (normalised) cumulated open mag- 
netic flux, integrated from the equator to a co-latitude 6. The 
curves indicate that 50% of the total open magnetic flux of a 
hemisphere is on average located below/above ~ 45°, with vari- 
ations of ±15° depending on the observational epoch. These 
results show that the distribution of open magnetic flux can be 
distinctively different from spot distributions determined from 
surface brightness maps alone. 

The present work extends the study o flSolankietal] ill 9971) 
and investigates the impact of latitude-dependent flux distri- 
butions on the rotational evolution of solar-type stars in more 
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Fig. 1. Cumulative open magnetic flux distributions, Og/Oo 
(cf. Sect. 12. 21 . on the visible hemisphere of AB Dor (P = 
0.5 1 d) , based on ZDI observations sec u red between 1995- 
2003 (bonati & Collier Cameronl 1 19971 bonati et alJ 1 19991 
120031) . Thin curves: 1995 (650% = 54°, solid); 1996 (47°, long 
dashed); 1997 (37°, short dashed); 1998 (49°, dashed dotted); 
1999 (60°, dotted). Thick curves: 2000 (29°, solid); 2001 (53°, 
long dashed); 2002 (56°, short dashed); 2003 (63°, dashed dot- 
ted). 



detail. Using the magnetic wind model described in lHolzwarthl 
l l2005l) . we quantify the influence of prescribed flux distribu- 
tions with different degrees of non-uniformity on the rotational 
evolution of stars, and verify their importance by comparing 
their impact with the influence of other magnetic field-related 
model parameters. 

2. Model setup 

The stellar structure of the low-mass star considered here con- 
sists of an outer convective envelope and an inner radiative 
core, each taken to be in solid-body rotation with possibly dif- 
ferent rotation rates. With J = IQ. being the AM, I the moment 
of inertia, and O the rotation rate, the rotational evolution of the 
star is determined by the set of coupled differential equations 



d\gt 
dig Q c 
dlgt 



digh 
dlgt 

d\gl c + _f 
dlgt J, 



t t ■ 



+ Je.Rl + J\ 



+ 1 



c.RIj 



(1) 

(2) 



which comprise changing moments of inertia and the AM 
transfer across appropriate boundaries; indices 'e' and 'c' de- 
note quantities of the envelope and of the core, respectively. 

2.1. Internal angular momentum redistribution 

The stellar structure is taken to be spherically symmetric. The 
temporal change of the moments of inertia is determined by an 
evolutionary sequence of stellar models of a 1 Mq star, which 
was generate d with an updated v ersion of the stellar evolution 
code of Kippenhah n et alJdl967l) . Figure|5]shows the evolution 
of the (outer) radii, masses, and moments of inertia of both the 
convective envelope and the radiative core. 

As the core-envelope interface recedes outward, the dy- 
namical stability properties at the base of the convection zone 
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Fig. 2. Evolution of the outer radii, r e / c (Panel a), masses, M e / C 
(Panel b), and moments of inertia, 7 e / c (Panel c) of the convec- 
tive envelope (solid) and radiative core (dashed) of a 1 M0 star. 
For t < 1 .5 Myr the star is fully convective. Its MS phase starts 
at t st 40 Myr; the vertical dashed line marks the age of the Sun. 



change and originally unstable mass settles down on the radia- 
tive core. This mass transfer is accompanied by an AM transfer, 



(3) 



The temporal change of the core mass, M c , is determined from 
the sequence of stellar models (Fig. 05); since the mass loss of 
the envelope through the stellar wind is negligibly small, it is 
M e = -M c . 

Hydrodynamic and hydromagnetic interaction at the 
boundary between the radiative and convective regions 
are expected to entail a coupling between the core and 
the envelope. Different coupling mechanisms based on 
magneto-viscous interaction or l arge-scale internal circula - 
tions have been investigate d by Pinsonneault et alJ 
Charb onneau & MacGreg orl dl993l) : iBouvier et alJ 
Allain (1998) and references therein. Here, a more quali- 
tative core-envelope cou pling model is used, f o llowin g the 
parametric approach of iMacGregor & Brennerl Jl99ll) and 
Kep pens et al.1 d 19951) . Based on the Rayleigh criterion a dy- 
namical stable state of rotation inside a star requires 1 the in- 
crease of the specific AM with increasing radius, d(£lr 2 )/dr > 
0. Given that the rotation rate of the convective envelope is 
braked and smaller than the rotation rate of the core, the sta- 
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in the equatorial plane and in the absence of viscosity 



4 



V Holzwarth and M Jardine: A further 'degree of freedom' in the rotational evolution of stars 



bility condition is violated at the core-envelope interface and a 
rotational instability sets in, transferring AM, 



Je.Rl - 



A./ 



-Jc. 



(4) 



from the core to the envelope to eliminate the differential rota- 
tion between the two regions. The AM required to equalise the 
two rotation rates, 



AJ 



L+L 



(a - a) , 



(5) 



is transferred on a timescale, t c , which is supposed to char- 
acterise the various visco-magnetic coupling mechanisms. The 
coupling time quantises the possibility to deposit AM in the 
core during the PMS phase and its retarded transfer to the en- 
velope in the course of the MS evolution. For the sake of sim- 
plicity the value of t c is taken to be constant during the entire 
rotational evolution of the star. 

2.2. Latitude-dependent magnetised stellar winds 

The AM loss rate of the convective envelope through a latitude- 
dependent magnetise d wind is determined following the ap- 
proach of lHolzwarthl J2005I) . The stationary, polytropic stellar 
wind is assumed to be symmetric with regard to both the rota- 
tion axis and the equator. The poloidal component of the mag- 
netic field is taken to be radial, with field lines forming spi- 
rals around the rotation axis on coni with constant opening an- 
gles. The whole stellar surface contributes to t he wind, wit hout 
'dead zones' retaining mass from escaping (cf. Mestel & Spruit 
119871) . 

The magnetic wind properties are determined through the 
radial magnetic field strength, 



B r , (t, 0,fle) = 



r (t ) 
ro(t) 



(fi<+Afi(Qe) •/(£)) 



(6) 



given at a reference level, ro, close to the stellar surface. The 
time-dependent radius ratio is to ensure that the total magnetic 
flux, 



In 7T 



<M«e) 



= J j B n0 r% sin ed6d$ = 47rrl(t )B r ,o (t , Q. e ) 



(7) 



o o 



only depends on the rotation rate of the star ISaarll 19911) . and 
an arbitrary reference time, to. The efficiency of the underlying 
dynamo mechanism is expected to increase with the rotation 
rate of the convection zone. The field strength variation, AB, is 
therefore determined in a way that the surface averaged field 
strength, 



B r ,o(to,Q e ) 



x/2 

= $B rfi 



sin Odd , 



(8) 



obeys the functional behaviour shown in Fig. [5] The non- 
uniform flux distributions, superposed on a constant 'back- 
ground' field, B K , are characterised by an enhancement of mag- 
netic flux at non-equatorial latitudes (Fig.@}: 
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Fig. 3. Dependence of the surface averaged radial magnetic 
field strength, B,-o(?o), on the rotation rate, Q e , of the convec- 
tive envelope. For large rotation rates the dependencies fol- 
low approximately linear (solid), sub-linear (long dashed), or 
super-linear (short dashed) power laws, oc (Q e /Q e0 )" n , with 
nn = 1,0.75, and 1.25, respectively. The deviation from the 
power law at small rotation rates is due to the constant back- 
ground field strength. 



AB 





Fig. 4. Latitude-dependent magnetic field distributions, f(0), 
close to the stellar surface. The non-uniformity of the Coronal 
Hole (CH) and Latitudinal Belt (LB) model is parametrised 
through the co-latitude 9q. Dashed lines show the lower, B K , 
and upper, B> = B< + AB, field strengths of each field distribu- 
tion. 



Coronal Hole (CH) model 

(I forO<0<0 o 
■' w \ for 6»o < 8 < 90° 

Latitudinal Belt (LB) model 

f(6) = cos 16 (9 - 6 ) 



(9) 



(10) 



Whereas the Latitudinal Belt model closely resembles what 
is found from field extrapolations based on ZDI images 
dMcIvoretalJ f2004). the Coronal Hole model is motivated 
by nu merous DI surface brightness maps (e.g., IStrassm eier 
120021) . Different degrees of high-latitude flux concentrations 
are realised by changing the co-latitude Bo of the analyti- 
cally prescribed functions / (Fig.|5J- We classify non-uniform 
flux distributions by the location of their 50%-open flux 
level, 650%, where the cumulated open magnetic flux, (Dg = 
^0^0 sm 6'dO', reaches half of the total value, Oq. 
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Fig. 5. Cumulative open magnetic flux, (S)g /<f>o, of the Coronal 
Hole (solid, for O = 45°) and Latitudinal Belt (dotted, for O = 
50°) model. Respective curves for a constant (long dashed) and 
dipolar (short dashed) field distribution are show for compari- 
son; see also Fig.[0for observational results. 

The total, surface integrated AM loss rate, 

/w^wo^J(^) 4 (|)(^)sin 3 ^, (11) 

o 

is expressed in terms of the plasma density, p A , and radial flow 
velocity, v r ,A, at the Alfvenic point, r A , wher e the flow velocity 
of the wind equals the Alfven velocity ("cf. lHorzwarthl l2005 ). 
Whereas p A , vv,a, and r A are functions of the co-latitude, 0, 
and subject to the latitude-dependent field distributions, Eq. 
ijSJl, the respective quantities, Pa,v k a< an d r A, are determined 
in the equatorial plane (0 = n/2) using the surface averaged 
field strength, Eq. iJS}. The latter quantities determi ne the AM 
loss rate following the approach of lWeber & Davisl i 19671) . 

8tt 4 

Jwd = —D.r^pAVr.h , (12) 

which is based on the simplifying assumption that the equato- 
rial wind structure can be generalised to all latitudes. 

2.3. Reference model parameters 

The wind structure is determined through boundary conditions 
prescribed at the reference level ro(f) = r e (t) + Ar, over the 
range of co-latitudes < < n/2, measured from the stellar 
north pole. r e (= R*) is the time-dependent outer radius of the 
convective envelope (Fig. |2^), and Ar = 0.1 Rq a constant ra- 
dial offset to locate the reference level of the wind at the base 
of the corona. 

For cool stars other than the Sun thermal wind properties 
are poorly constrained by observations; for possible constraints 
resulting from relationships between the temperature and den- 
sity of closed coronal loops and rot ation/Rossb y number see, 
for ex ample, Llordan & Montesinosl dl99ll) : llvanova & TaarrJ 
J2003I) . In the following solar-like values are assumed for the 
temperature, Tq = 2 ■ 10 6 K, and (particle) density, «o = 
10 8 cm 4 . The entropy change of the wind with increasing dis- 
tance from the star is quantified through the polytropic index, 

r = i.i5. 



Table 1. Mean deviation, (Q c - Q e ) /Q c , from isorotation at so- 
lar age, t Q = 4.7 Gyr, for different coupling timescales, t c . 

r c [Myr] 1 5 10 15 20 50 100 
Afi/n c [%] 0.3 1.3 2.6 4.1 5.6 16.4 35.7 



The surface averaged magnetic field strength, defined in 
Eq. l|8}, is taken to follow the quasi-linear power law shown 
in Fig. (solid line), that is B rj0 = (1 + 1.5 (O e /f2 )) G, with 
Q.q = 2.8 • 10~ 6 s _1 . This is in agreement with observations of 
rapidly rotating stars, which show magnetic field strengths up 
to ab out two orders of m agnitude larger than in the case of the 
Sun JDonati & Collier CamerorJl997h . 

Helioseismological observations sh ow that the rotation rat e 
in the solar interior is roughly uniform JThompson et alJ2003t) . 
The rotation rate of the solar radiative core and convection zone 
are within ~ 4% about £2 = 2.8 • 10 _6 s _1 . The present Sun 
thus constrains the coupling timescale to the effect that the 
value of t c ought to achieve isorotation within a few percent 
at the solar age, f© ~ 4.7 Gyr. Using the reference model pa- 
rameters described above, we accomplished simulations with 
different coupling timescales to determine the relative devia- 
tion, (Q c - Q e )/Q c , from isorotation at solar age (Table 
For t c , ^ 15 Myr, a value si milar to the one ad opted by 
iMacGregor & BrenneJJl99ll) orlkeppens et alJ l l 1995b . the de- 
viation is in accord with the observational constrains given by 
the Sun. 

3. Results 

We determine the rotational evolution of 1 Mq stars with initial 
rotation rates SIq = 5 ■ 10~ 6 , 2 ■ 1CT 5 , and 6 • 1CT 5 s _1 from the age 
1.6 Myr onward. The rotation periods, between 1.2 d and 15 d, 
approximately span the observe d range of rotation p eriods of 
young (< 5 Myr) T Tauri stars jBouvier et alJll993i) . In their 
initial state the convective and radiative zones are taken to be 
in isorotation, that is Q c = fi e = Oq. 

The rotational histories determined with the AM loss rate 
/w = -Avd. that is using the WD approach, Eq. H2\ . are used 
as reference cases (Fig.|6j. In the course of the PMS evolution 
the rotation rate of the envelope is dominated by its decreasing 
moment of inertia and the AM loss which goes with the mass 
settling down onto the radiative core, whereas during the MS 
phase only the magnetic braking and the internal coupling are 
relevant. 

3.1. Non-uniform flux distributions 

We determine the rotational history of stars subject to the 
Latitudinal Belt (LB) and Coronal Hole (CH) flux distribu- 
tions with 50%-open flux levels located at different latitudes. 
At first, the non-uniform flux patterns are taken to be station- 
ary (#50% = const.) during the entire evolution, to separate their 
influence from other rotation-dependent effects. 

An accumulation of open magnetic flux at high (low) lati- 
tudes causes, with respect to an uniform flux distribution, a re- 
duction (enhancement) of the AM loss rate through the stellar 
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Fig. 7. Relative deviations of the envelope rotation rate for the Latitudinal Belt (left) and Coronal Hole (right) flux distributions. 
For 50%-open flux levels at low latitudes (i.e., large co-latitudes 650%, labels) the deviations are negative (dashed), and for high 
latitudes positive (solid). 
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Fig. 6. Rotational evolution of the convective (thick) and radia- 
tive (thin) zone of stars subject to an AM loss rate according to 
Web er & Davisl (119671) . The dashed line indicates the break-up 
rotation rate, for which the co-rotation radius equals the actual 
stellar radius, fibu = ^GM^/rl. The cross marks the rotation 
rate of the present Sun. 

wind. At the age of the present Sun the resulting deviations of 
the stellar rotation rate are found to cover a range of values be- 
tween about -40% and 200% different than the respective ref- 
erence cases (Fig.0. The relative deviations show two distinc- 
tive regimes, corresponding to the PMS and the MS phase. For 
the deviations to increase (decrease) with time the net AM loss 
rate of the envelope has to be smaller (larger) than in the ref- 
erence case. Since we consider the difference between two ro- 
tational histories, the rotation-independent mechanisms cancel 
out, so that only the magnetic braking and the internal coupling 
are relevant. The magnetic braking is a direct consequence of 
the present stellar rotation and causes an immediate AM loss of 
the envelope. The internal coupling, in contrast, is proportional 
to the differential rotation and during the early evolution rather 
inefficient, due to the initial isorotation (Fig. [8}. In the PMS 
phase the difference between the AM loss of the envelope and 
its internal AM gain is consequently large and the deviations of 
the rotation rate quickly increasing. After arrival on the MS the 
differential rotation is adjusted to the new situation and the tap- 
ping of the AM reservoir of the more rapidly rotating core now 
replenishes most of the AM loss carried away by the magne- 
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Fig. 8. Relative timescales (with respect to respective reference 
cases) on which the rotation rate of the envelope changes due 
to the magnetic braking (solid) and internal coupling (dashed), 
for Qrj = 2 • 10~ 5 s and a LB distribution with f5go% = 30 and 
75°. 

tised wind; the further increase of the deviations is respectively 
weaker. 

The actual spin-down time scale of the envelope is longer 
than the braking time scale, because the AM gain of the en- 
velope through the internal coupling follows closely its AM 
loss through the magnetic braking (Fig. [8}. Note that in the 
course of the evolution the latter is found to converge to- 
ward the reference value. Since the magnetic field strength de- 
creases with the rotation rate the relative contribution of the 
(latitude-dependent) magneto-centrifugal driving to the overall 
wind acceleration becomes smaller and the one of the (latitude- 
ind ependent) th ermal driving larger. The difference between 
the lWeber & D avis and the present non-uniform wind approach 
thus becomes smaller and the spin-down timescale similar. 

The non-uniform flux distributions with f5 50 % = 60° have 
nominally the same 50%-fiux level as the uniform flux distri- 
bution of the reference cases. The deviations of the relative ro- 
tation rate of ~ 10% in the case of the LB model indicate that 
the usage of 6*50% only allows for a rudimentary classification 
of flux patterns. By comparing the rotational histories resulting 
from different flux distributions with equivalent 50%-open flux 
levels, we find that over the range 6> 50 % = 10 - 60° the peaked 
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Fig. 9. Relative deviations {dashed lines indicate negative val- Fig. 10. Relative deviations in the case of non-linear dynamo 
ues) in the case of rotation-dependent 50%-open flux levels efficiencies, AB oc Q"*. For = 0.7,0.8, and 0.9 (solid, top 
#50% = 90° (Q/Q0)"" 1 , with rini = -0.25, -0.5, and -1 {labels), down) the deviations are positive, and for «d e = 1-3, 1.2, and 



Its value is constrained to the range 15° < #50% < 60° 



1.1 (dashed, top down) negative. 



Latitudinal Belt and the bi-modal Coronal Hole model yield 
rotational evolutions which are within a tolerance of 15% con- 
sistent. In the special case of a dipolar field distribution (with 
f(ff) = cos 9 and f? 5 o% = 45°, Fig. [5}, the correspondency is 
found to be good within about 10%. 

Observations show that surface distributions of starspots 
depend on the stellar rotation rate, with spots being prefer- 
entially located at higher latitudes the faster the star rotates. 
We investigate this aspect by assuming a rotation-dependent 
latitude of the 50%-flux level, which follows the power law 
5O% = 90°(n/a Q ) n " 1 , with n rd = -0.25,-0.5, and -1; this 
simple relationship is only used to examine the basic effects 
and not meant to reproduce any complex observed or theoret- 
ically derived relations. To avoid unrealistic high flux concen- 
trations in the vicinity of the poles or the equator, the 50%- 
flux latit udes are const rained to the range 15° < #50% < 60° 
(cf. lSolankietalJl997l) . Figure|9]shows the relative deviations 
of the rotation rates from the reference values in the case of a 
rotation-dependent LB distribution; for the CH model the re- 
sults are very similar. The major consequence of the rotation- 
dependence is the enhanced stellar spin-up and spin-down dur- 
ing the PMS and the late MS phase, respectively. The evolu- 
tionary stage at which the rotational history is altered by the 
shifting of magnetic flux to lower latitudes depends however 
on the actual functional dependence of # 50% (Q). 

3.2. Comparison with dynamo efficiency and 
saturation 

To examine whether the impact of non-uniform flux distri- 
butions is significant with respect to other magnetic field- 
related model assumptions, we compare the deviations de- 
scribed above with those obtained by successively changing 
the functional dependence of the magnetic field strength on the 
stellar rotation rate in the framework of the reference cases (i.e., 
with uniform flux distributions). 

A linear dynamo efficiency is an acceptable approximation 
in the case of slowly rotating stars, but it is found to fail for 
fast rotators. We analyse the sensitivity of the rotational evolu- 
tion on changes of the dynamo efficiency by re-calculating the 



reference cases, but now subject to non-linear magnetic field- 
rotation relations, AB oc Q" dL , with «de + 1 (cf. Fig. [3}- During 
the MS phase, the resulting relative deviations cover a large 
range between -50% and 150% (Fig. 1 1 0i. Sub-linear (super- 
linear) dynamo efficiencies imply a weaker (stronger) increase 
of the magnetic field strength with rotation rate and conse- 
quently a moderation (enhancement) of the magnetic braking. 
Their influence is therefore qualitatively similar to rotation- 
dependent concentrations of magnetic flux at high and low lat- 
itude, respectively. To determine explicit relations between the 
dynamo efficiency and the location of the 50%-open flux level, 
we generate grids of rotational histories with different «d e - and 
#50% -values, and associate corresponding values by comparing 
the relative derivations of the rotation rate from the reference 
case 2 . Figure[^shows quadratic fits of the resulting grid-based 
relations for the LB and CH flux distribution. The relationships 
show that non-uniform flux distributions can imitate a large 
range (here: between -45% and 15%) of dynamo efficiencies 
and, consequently, if not taken properly into account, conceal 
them from an observational determination. 

The limiting field strength beyond which dynamo satura- 
tion occurs is referred to in terms of a critical rotation rate, Q sat , 
that is AB(Q. e > Q sat ) = AB(£2 sat ). We repeat our calculations of 
the reference cases under the assumption of saturation rotation 
rates in the range of £2 sat = 2 - 80 Q (see Fig.[21for exam- 
ples). The saturated magnetic braking causes an enhanced spin- 
up during the PMS phase, higher rotation rates on the ZAMS, 
and a weaker spin-down during the MS phase until the rota- 
tion rate descends into the regime of the non-saturated dynamo. 
We compare the rotational histories resulting from a saturated 
magnetic braking with those resulting from non-uniform flux 
distributions to associate 50%-open flux level with correspond- 
ing dynamo saturation rates (Fig. 1131 . For given initial stellar 
rotation rates and flux distributions, the curves 3 describe nomi- 
nal dynamo saturations which cause similar rotational histories 

2 Since the variations of the rotational histories over the entire evo- 
lutionary time range are qualitatively different, we constrain the com- 
parison to the early MS phase between 50 Myr <t< 800 Myr. 

3 Note that this comparison considers exclusive mechanisms only, 
that is how a saturated dynamo with a uniform flux distribution cor- 
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Fig. 11. Relationships between the dynamo efficiency, n^, and 
the latitude, 650%, of the 50%-open flux level in the case of 
the Latitudinal Belt (LB, top) and Coronal Hole (CH, bottom) 
distributions. The initial stellar rotation rates are Qo = 6 ■ 
(solid, crosses), 2- 10 5 (long dashed, triangles), and 5 • 10~ 6 s -1 
(short dashed, rhombs). 
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Fig. 12. Relative deviations of the envelope rotation rate in 
the case of a dynamo saturation beyond Q sat = 10 (solid), 
20 (long dashed), 30 (short dashed), 40 (dashed-dotted), and 
50 O (dotted). 



like a non-uniform flux distribution with a 50%-open flux level 
at the respective co-latitude. If the actual dynamo saturation 
rate is larger than then nominal one, the rotational evolution 
is likely to be dominated by the impact of the non-uniform 
flux distribution. If, in contrast, the actual dynamo saturation 
rate is smaller than the nominal one, then the behaviour of 
the magnetic dynamo is expected to dominate the evolution. 
A concentration of magnetic flux at high latitudes of rapidly 
rotating stars can therefore imitate a rotational evolution like 
a dynamo saturation limit of about 30 - 40 fi - In case a con- 
siderable amount of open flux is located at intermediate or low 
latitudes, as indicated by recent field reconstructions based on 
ZDI observations, the nominal saturation level would be rather 
high (> 60 Q.q) and therefore probably beyond the actual dy- 



responds to an unconstrained dynamo with non-uniform flux distribu- 
tion. 



Fig. 13. Relationships between the nominal saturation rotation 
rate, £2 sat , and the latitude, #50%, of the 50%-open flux level for 
stars with initial rotation rates Qo — 6 ■ 10~ 5 (solid), 2 ■ 10~ 5 
(long dashed), and 5 • 10~ 6 s~' . In contrast to the Coronal Hole 
model (thin lines), Latitudinal Belt distributions (thick lines) 
with 050% > 53° have no corresponding dynamo saturation. The 
shaded stripes indicate the rotational regimes where the chro- 
mospheric and coronal emission (dark grey) and the variation 
of the photometric light curve (light grey) are observationally 
found to saturate. 

namo saturation limit. In this case non-uniform flux distribu- 
tions alone could not account for the very high rotation rates of 
young stars. 

4. Discussion 

Our results show that a concentration of open magnetic flux at 
high latitudes of active stars yields rotational histories which 
deviate up to about 100% (early MS phase) to 200% (late 
MS phase) from the rotation rates obtained in the case of uni- 
form surface fields. The reduction of the AM loss rate entails 
a quicker increase of the stellar rotation rate during the PMS 
spin-up and a moderated spin-down during the MS evolution. 
The higher the (initial) rotation rate, the larger the AM loss 
due to magnetised winds, and the larger the susceptibility to 
variations in the open flux distribution. The influence on the 
rotational history is a cumulative effect, which depends on the 
evolutionary stage of the star but eventually also on the applied 
model assumptions. 

4.1. Model considerations 

In our wind model the poloidal magn etic fi eld component is 
prescribed to be radial (cf. lHolzwarthll2005h . whereas a fully 
consistent treatment of the multi-dimensional problem, includ- 
ing the trans-field component of the equation of motion, is 
found to show a collimation of open field lines towar d the ro- 
tation axis with increasing distance from the star (e.g. JSakurail 
119851) . The influence of this effect on our results is difficult to 
quantify, because investigations focused on this phenomenon 
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are usually limited to particular and illustrative cases. Another 
approach to include a n on-uniform magnetic fie ld topology is 
the stellar wind model o f lMestel&Spruitl i ll 9871) . which incor- 
porates in the vicinity of the stellar surface polar 'wind' and 
equatorial 'dead zones', the latter preventing plasma escaping 
from the star and thus reducing both the mass and AM loss 
rate. Whereas in this model the extent of the dead zone depends 
on the rotation rate of the star, the magnetic field distribution 
at the sur face is prescri bed to be dipolar. Based o n the for- 
malism of JMeste]l fl968) and lMestel & Spruilllfl987i) . lKawaie3 
Jl988l) derived a parametrised description of stellar AM loss 
rates through stellar winds, which also comprises qualitative 
variations of the magnetic field topology; a definite association 
of non-uniform magnetic flux distributions with the respective 
model parameter is however missing. Although it remains to 
be investigated how more complex flux distributions alter the 
magnetic field topology of stellar winds including effects like 
collimation or dead zones, we deem these aspects to be less 
crucial since our analyses are based on relative deviations be- 
tween two rotational histories. However, more detailed studies 
concerning this hypothesis are required. 

The surface averaged magnetic field strengths, cover- 
ing here roughly two orders of magnitude over the range 
of relevant rotation rates, are consistent with observations 
both in the case of the Sun and rapidly rotating stars (e.g., 
iDonati & Collier Cameronl Il997t) . The localised peak field 
strengths are, depending on the location of the 50%-open flux 
level and the non-uniformity of the flux distribution, even larger 
(up to the order of kilo-Gauss), which is also in agreement with 
observations. 

The even flux distribution of the Coronal Hole model and 
the peaked distribution of the Latitudinal Belt model simu- 
late complementary non-uniformities and allow us to assess the 
sensitivity of the results beyond our 50%-open flux classifica- 
tion. Whereas both cases are consistent within < 15%, we ex- 
pect this tolerance level to be larger for more complex flux dis- 
tributions. But in view of the current observational and model 
limitations a more sophisticated classification scheme appears 
yet to be inappropriate. 

The thermal wind properties are described by solar-like val- 
ues, which are taken to be independent of the rotation rate, the 
stellar latitude, and the evolutionary stage of the star. Although 
there are indications for a rotation-dependenc e of the coro- 
nal temperature and de nsity (e.g. Jjordan & Montesinoslfl99lt 
llvanova & Taaml Eo03). it is difficult to derive accurate wind 
parameters, since observational wind signatures are veiled by 
the outshining coronal emission. The polytropic index of T = 
1.15 implies an efficient heating and thermal driving of the stel- 
lar w ind; its value lies well in between those of s imilar studies 
(e.g.. lSakurailll985UKeppens & GoedbloedboOOl) . It has been 
chosen to ensure stationary wind solutions at all latitudes even 
for small rotation rates, when the magneto-centrifugal driving 
is inherently small. This implies that the contribution of the 
latitude-independent thermal driving is relatively large com- 
pared with the latitude-dependent magneto-centrifugal driving. 
In case of a weaker thermal driving, either because of cooler 
coronae or a smaller energy flux in the wind, the impact of 
non-uniform flux distributions is expected to be even stronger 



than described above. If, in turn, magnetised winds turn out to 
be hotter, latitude-dependent effects may be diminished (given 
that the thermal wind properties themselves are not latitude- 
dependent). 

The rapidly rotating stellar core represents an AM reservoir 
which is tapped by the envelope during the MS evolution. The 
efficiency of the internal coupling, quantified through a char- 
acteristic timescale (here determined to be r c = 15Myr), is 
under debate, with previous studies investigating values over 
a very large r ange, from very short (solid-body rotation; e.g., 
iBouvier et alJll997|) to very long (essentially decoupled differ- 
ential rotation: Ijianke & Collier Cameronlfl993t lAllainll 1998b 
timescales. The results are typically not fully consistent with 
all observational constraints given by the rotational distribu- 
tions of stars in young open clusters of different age, which 
may point out possible deficiencies of previous models. In this 
respect the non-uniformity of surface magnetic fields presents 
an additional and so far rather ignored 'degree of freedom', 
which may help to match theoretical and observed rotational 
histories of stars. 

4.2. Relevance to dynamo efficiency and saturation 

Our comparison of the impact of non-uniform flux distributions 
with the influence of other magnetic model parameters under- 
lines its importance. The dynamo efficiency, that is the depen- 
dence of the (here, open) flux generation on stellar rotation, is 
generally accepted to increase with the rotation rate. Analytical 
relations based on the theory of aQ-dynamos (e .g., between 
the dynamo number and the Rossby number; iNoves et alJ 
Il984t iMontesinos et alJEool as well as empirical relation- 
ships between the rotation rate and magneti c activity s ignatures 
like the coronal EUV and X-ray emiss ion ( Mathioudakis et alJ 
Il995l iHempelmann et alJ Il995t IPizzolato et al.1 120031) imply 
power law- dependencies over a large range of rotation rates. 
Skumanich ( 1972) found the evolution of chromospheric ac- 
tivity to decrease with time following Q oc r 1 ' 2 . This rela- 
tion implies a linear increase of the characteristic stellar mag- 
netic field strength with rotation rate; an example for the com- 
monly agreed concept that the decrease of magnetic activity 
with age reflects the e voluti o n of the stellar AM and rota- 
tion rate. According to ISaarl J 1 99 lit , it is however more the 
magnetic flux which increases linearly, instead of the field 
strength. The linear dynamo efficiency is widel y used in stud- 
ies about the rotational evolutio n of stars (e.g.. lKeppens et all 
h995tlKrishnarnirrthi et alll~997l) . Here, for deviations from the 
linearity between ±30% the rotation rates on the early MS are 
found to differ between 150% and -50%. Our analysis has 
shown that non-uniform flux distributions can easily produce 
similar deviations. A large fraction of open flux in the vicinity 
of the stellar poles can therefore compensate the influence of 
a super-linear dynamo efficiency, feigning a less efficient (e.g., 
linear) dependence on the rotation rate. 

A concentration of magnetic flux at high latitudes has qual- 
itatively a similar influence on the rotational history as a dy- 
namo saturation. Since both effects entail a reduction of the 
AM loss rate, their synergistic action is expected to enable even 
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higher rotation rates. Figure [O] allows for estimates (based 
on the comparison of exclusive effects) of the principal con- 
tribution to the overall reduction of the AM loss. If the sat- 
uration limit of the chromospheric and/or coronal emission 
(F sat uv/x ~ 3 ... 1.5 d) is representative of dynamo saturation, 
then we expect non-uniform flux distributions to contribute 
only marginally to the reduction of the AM loss. If, however, 
the photometric saturation limit (-P sat ,phot 5 12h) reflects the 
saturation of the dynamo mechanism, then the rotational evo- 
lution of most stars is dominated by the actual non-uniformity 
of the surface flux distribution. But the saturation of indirect 
activity signatures does not a priori indicate a saturation of the 
underlying dynamo processes in the convective envelope, since 
the former also depend on, for example, atmospheric proper- 
ties and radiatio n processes which follo w different functional 
beha viours (e.g..|Unruh & Jardinel 19971) . 

ISolankiet alJ i 19971) applied a bi-modal flux distribution, 
characterised through polar flux concentrations, to the rota- 
tional evolution of stars. For somewhat different model param- 
eters than ours they found a resemblance between the rotational 
histories subject to their coronal hole model and a dynamo satu- 
ration limit of 20 Q , respectively. Since this saturation limit is 
of the order required to explain the presence of rapidly rotating 
stars (~ 10 - 20 Q ), they question the necessity of a dynamo 
saturation at low rotation rates and argue in favour for values 
> 50 Q . Whereas our investigation confirms their results in 
principle, we find that reductions of the AM loss rate resulting 
from high-latitude flux concentrations are smaller; according 
to our model a flux distribution with, for example, 850% * 15° 
is equivalent to a dynamo satu ration limit o f about 40 £2 . The 
difference with respect to the ISolanki et alJ value (~ 2O£2 ) 
may be due to different model assumptions and parameters. 
Since our value is still below their supposed dynamo saturation 
limit of ~ 50 £2 , a strict concentration of magnetic flux around 
the poles would be able to dominate the rotational evolution 
of stars. Observations indicate, however, that the 50%-open 
flux level is likely located at somewhat lower latitudes. Surface 
brightness maps occasionally show and dark elon gated features 
reaching down to intermediate latitudes ( Strassmeiei 2002, and 
references therein), which shift the mean 50%-flux level (av- 
eraged over longitude and evolutionary timescales) equator- 
wards. But the presence of magnetic flux in the form of dark 
spots is a priori not equivalent with the presence of open mag- 
netic field lines. A more striking constraint arises from ZDI ob- 
servations, which in combination with field extrapolation tech- 
niques allow for the reconstruction of the magnetic field topol- 
ogy and the determination of the actual surface and latitudinal 
distribution of open magnetic fields. Although the number of 
detailed observations is yet rather small, first results show that 
the principal part of (dete ctable) open flux is apparently located 
at intermediate latitudes (Mclvor et al. 2004), which motivated 
our Latitudinal Belt model. The cumulated flux distributions of 
AB Dor (cf. Fig.Q show that the average 50%-flux level is in- 
deed located between 30° < 9^% < 60°. In this range of values 
our equivalent saturation limit for rapid rotators is > 60 Q and 
the impact of the non-uniform flux distribution consequently 
smaller than the influence o f the d ynamo saturation at ~ 50 £2 
supposed bv ISolanki et alJ 09971) . In this sense, we consider 



the effect of non-uniform flux distributions more as a comple- 
mentary rather than an alternative mechanism for the formation 
of rapid rotators, which is however based on observationally 
verifiable principles. 

4.3. Observational constraints 

Compared to the large number of obs erved surface brigh tness 
distributions of rapidly rotating stars ( Stras smeien l2002). ZDI 
observations of actual magnetic flux distributions are yet rather 
sparse. A larger database in this field is required to put tighter 
constraints on possible open flux distributions, that is the quali- 
tative and quantitative description of its non-uniformity and de- 
pendence on the stellar rotation rate and the evolutionary phase 
of the star (including its dependence on stellar mass). Whereas, 
for example, brightness distributions indicate a clear poleward 
displacement of starspots with increasing rotation rate, for open 
flux a similar relationship is yet not available. ZDI observa- 
tions are handicapped through insufficient signal levels from 
dark surface regions, so that an unambiguous determination 
of the flux topology inside dark polar caps (i.e., whether it is 
mainly unipolar and open or more multi-polar and clos ed) is yet 
hardly possible JPonati & BrowrJ 19971: [Mclvor et alJ2003l) . If 
the field reconstructions in the case of AB Dor and LQ Hya 
prove to be representative and characteristic for rapid rotators, 
then we consider the contributions from high-latitude regions 
to the total AM loss rate to be of minor importance for the ro- 
tational evolution of rapid rotators. 

Owing to their qualitatively similar behaviour it is question- 
able whether the impact of non-uniform flux distributions can 
be observationally separated from the influence of a dynamo 
saturation. This potentially limits the usefulness of rotation rate 
distributions as tests of dynamo theories, unless they are sup- 
plemented by more detailed information about the structure of 
stellar magnetic fields. A larger observational database is re- 
quired to refine the rotation rate beyond which the location of 
the 50%-open flux level is high enough to cause a discernible 
impact on rotational distributions of cluster stars. Supposed that 
the limiting rotation rates for high-latitude 50%-flux levels and 
dynamo saturation are sufficiently low and high, respectively, 
it may be possible to distinguish characteristic signatures of 
non-uniform flux distributions in the (differential) distribution 
of intermediate rotators with rotation rates in the domain con- 
strained by the two limiting values. 

An important role may fall to theoretical models concern- 
ing the pre-emptive evolution and post-eruptive transport of 
mag netic flux to high latitudes in rapidly ro tating young stars 
(e.g. JGranzer et alJ2000l:lMackav et alJ2 004). since they allow 
for a verification of our insight into stellar magnetic proper- 
ties by means of activity proxies like (the rotational modulation 
of) co ronal X-ray and chromospheric UV emission JVilhu et alJ 
Il993t) . Observations in the UV/EUV spectral range may also be 
used to refine our assumptions about the thermal wind proper- 
ties, comprising its latitudinal and rotational dependence, be- 
cause the total AM strongly relies on the stellar mass loss rate 
and on the acceleration of the wind through thermal driving. 
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5. Conclusion 

Non-uniform magnetic flux distributions have a significant im- 
pact on the rotational evolution of stars, mimicking the effect 
of a large range of dynamo efficiencies and saturation lim- 
its. They present an additional degree of freedom in the mod- 
elling of stellar rotational histories, which can generate dif- 
ferences cumulating up to 200%. Neglecting their effect im- 
plies considerable uncertainties in other magnetic field-related 
model parameters such as the dynamo efficiency of up to about 
40%. Althoug h our results a re, in principle, in agreement with 
those o fSola nkiet all ill 9971) . we find the effect of non-uniform 
flux distributions to be less efficient than in their investigation; 
whereas they find that a concentration of magnetic flux at polar 
latitudes imitates a nominal dynamo saturation limit at 2OQ , 
we find values of about 35£2 . 

Non-uniformities in the form of strong flux concentrations 
at high latitudes efficiently reduce the AM loss through magne- 
tised winds, entailing high stellar rotation rates. Since magnetic 
field distributions reconstructed on the basis of ZDI observa- 
tions however indicate a considerable amount of open flux at 
intermediate latitudes, the anticipated reduction of the AM loss 
is expected to be smaller than implied by frequent DI observa- 
tions of high-latitude starspots and polar caps. The influence of 
non-uniform flux distributions alone thus appears to be insuf- 
ficient to explain the existence of very rapid rotators, but their 
significant moderation of the AM loss rate makes the require- 
ments for a dynamo saturation less stringent, enabling satura- 
tion limits > 40 Q . 
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